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' Abstract 

\ We consider a two particle system on a star graph with 5-function 

00 ' interaction. A complete description of the eigensolutions with real 

momenta is given; specifically it is shown that all eigensolutions can 
be written as integrals in the momentum plane of sums of products of 
1 appropriate one particle solutions. 

1 Introduction and Summary 

This paper is a continuation of [2] in which was considered the two par- 
| ticle problem with <5-function interaction on a star graph with n edges, 

' T n . We briefly recall the problem: the configuration space is = T n x 

r n with the local coordinate neighbourhoods consisting of n 2 quadrants 
{Qij = {(xi,yj)}}^j =l . The hamiltonian is defined as 



*2 a2 



with c real, A \ n = — ■£ — , and Kirchhoff boundary conditions at 

the boundaries of the quadrants 



^Qij^O - ^\Q kj ,x k =Q '■> ^2 Q Xl 



1=1 yl 



= a,Vi,j,k (2) 

Qij,x t =0 

= 0,Vi,j,k. (3) 

iu>vi=Q 



The 5 interaction on the diagonal D = {xi = yi ; Xi, in G Qii, i G {1,..., n}} 
means that eigensolutions satisfy the boundary conditions pQ 



1 / d± _ 8± 

2 \ dxi dyi 
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This problem (TjQ) is interesting because the equivalent problem on the line 
(n = 2) is 'completely solvable', not just for two particles but for an arbitrary 
number of particles [3J. In this paper we use the results from [2] and describe 
all solutions to the eigenvalue problem for H c with real momenta. To be 
precise we describe all eigensolutions of H c of energy A > which can be 
written in the form 

iP(x,y) = [ ( - ^x+iVx=Wy ^++ + e -ikx+iVx^y 
JVxi ^ 

- e -ikx-iV^k?y e ikx-iV>^y j dk (g) 

on the off diagonal quadrants, or on one of the sectors {x^ > in}, {xi < yi} 
of the diagonal quadrants. Here I = [0, 1] and since 



lim i f R f R \^\ 2 dxdy = 2^ [ |V 

u , / 



arl2 dk 



we assume that the transforms ip ar (k) are square integrable on yXI. We 
use the notation ipfj' to denote the transform on the quadrant Qy. 
In this paper we show that these (real momentum) eigensolutions of H c 
call all be written as integrals in the momentum plane of the 2n 2 — 2n 
eigensolutions of H c described in §4 of [2]. These solutions were characterised 
by the fact that they are finite sums of products of one particle solutions 
on the configuration space cut along the diagonal, -F^. Let us denote here 
this 2n 2 — 2n-dimensional vector space of solutions by W and the individual 
solutions by 

G l (h,k 2 ). 

In order to prove this we first consider solutions ([6]) satisfying conditions 
(0 [3j at the boundaries of the quadrants but which are not smooth on the 
diagonal and which did not, necessarily, satisfy (JHEJ). For want of a better 
term we refer to such solutions as 'basic solutions'. It is shown in $2] that 
(0 [3J imposes linear conditions on appropriate vectors of transforms i\) aT . 
Using these linear conditions we see that all 'basic solutions' can be written 
as integrals (in the momentum plane) of the n 2 + n + 1 solutions described 
in §3 of P]. In order to be precise and introduce some notation let us denote 
the solutions described in §3 of [2] by 

F i (k 1 ,k 2 ) 

and the n 2 + n + 1-dimensional vector space they span by V. Then it is 
shown in £}2]that the basic solutions correspond to L 2 (^f\I^ V under the 
map 

V / fiWF* (k,VX-kA dk (7) 
with (ft) G L 2 (v 7 !/) V. 

In ^3] we consider the imposition of the boundary conditions (JU E]) on the 
basic solutions. These impose linear conditions on the transforms, but only 
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if we consider the permutation of the one particle momenta, ie. we need 
to consider transforms at momenta k and Vl — k 2 , see (j 1 9 p : or equiva- 



lently F l yk, a/A — k 2 J and F l ^a/A — k 2 , k) ; or equivalently F 1 (ki, £2) and 

F l (k2,k\). On the other hand, in §4 of [2] all eigensolutions of H c (con- 
structed as products of one particle solutions) are found by enumerating 
those linear combinations of F l (fci,/^) an d F l (k2,k\) which satisfy the 
boundary conditions ([HE]). But this is exactly the problem described above 
from ^31 In this way we see that the set of all solutions, at given energy 
A > and with real momenta, to the eigenvalue problem for H c is simply 

the tensor product L 2 (^\J\l^j 8> W. The map from the tensor product to 

the eigensolutions of H c is 

Y, In; 9i(k) G l (k, a/A - k 2 ) dk (8) 



with ( gi ) e L 2 ( y I J J ® W. 

In the last part of 53 we recall the form of the 2n 2 — 2n eigensolutions G % . 
We present them here in a slightly more geometric 'nicer' form than in [2]. 

2 Description of the Basic Solutions 

Since the boundary conditions (0 [3]) are scale invariant we may rescale X{ 
and yi on T 2 thereby normalising A = 1 — in the following sections we put 
A = 1 in (0 El El). From © 

4>{x, 0) = / (e ikx [ - ^ ++ + + e~ %kx [</>_+ - V-]) dA; (9) 



(10) 



and 

so that the Kirchhhoff boundary conditions ([3]) are written 

for all i,j,k. Similarly, the Kirchhhoff boundary conditions ([2]) take the 
form 

for all i,j,k. These can be written in compact form by introducing £ 
M n (C 2 ) defined as 

£)•*-($)■ 
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Then (J31 [2]) can be written as 

& = - E > & = - E s * t **j • ( 12 ) 

fc Ac 

where S*j is the unitary matrix 2P—I with P the projection onto (1,1,..., 1)' 
and t is the transposition on C 2 . 

We denote by ip?j the transforms which, on the diagonal quadrants, give 
the eigenfunction in the sector X{ > yi. Similarly we denote by ipf^ the 
transforms which, on the diagonal quadrants, give the eigenfunction in the 
sector Xi < yf, in particular ijjfj' = ip°J for i ^ j. We define £, % an d £, 
X similarly. Then using (|12p we see that all solutions to the two particle 
problem satisfying (0 E]) (but not necessarily dHO) can be described by the 
system 

tt 1 ( X S - Stx) = , n 1 ( X - X) = (13) 

where ir is the projection onto the diagonal in M n . This equation is finite 
dimensional, and linear, implying that all basic solutions can be described 
by a finite dimensional vector space (tensored by a space of functions). 
It is shown in the appendix that the solutions of (|13[) can be split into four 
subspaces of dimension (n — l) 2 + 1, 2(n — 1), n — 1 and 2 respectively. 
Likewise the vector space V of solutions, from §3 of [2], splits into four: 
the [n — l) 2 + 1-dimensional smooth, symmetric subbasis; the 2(n — 1)- 
dimensional smooth, antisymmetric subbasis; the n — 1-dimensional non 
smooth, symmetric subbasis; and the 2-dimensional non smooth, antisym- 
metric subbasis. It is a straightforward, but long, calculation to show that 
under these four subbases of V are mapped onto the solution subspaces 
of (fT3j) of the same dimension. Consequently, the space of basic solutions is 
spanned by solutions of the form ([7]). 

3 Eigensolutions of H c 

The eigensolutions of H c may be defined as the subspace of basic solutions 
which satisfy the boundary conditions H |5| on D. In ^3.11 we continue 
the above analysis describing the eigensolutions in terms of the transforms 
ip°" r , ip (7T and showing that they are given by ([8]) . In ^3.21 we rewrite the 
eigensolutions described in [2] (the G l ) presenting them in a slightly more 
geometric way then in [2]. 

3.1 Boundary conditions on transforms on the diagonal 

The main point of this subsection is that, in analogy with (|13p . the boundary 
conditions (jUE]) may be described by a linear system on the transforms. To 
do this, however, we need to consider the transforms not just at momentum 
k but also at momentum Vl — k 2 , as in (1190 . This is equivalent to the fact 
that, in §4 of [2], we need to consider momenta (fci,^) and (k2,k\) for the 
first and second particles respectively. In this way we see that the linear 
system induced on (|19p by (JU [5]) is nothing more or less than the linear 
conditions on the diagonal satisfied by the solutions G 1 described in §4 of 
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On each diagonal quadrant we have from (|6|) 

il>(x = y + ) = [ ( - e Kk+VT^)x ^++ + e ^(k-VT^)x 



-i{k+VT^)x ^— + e ^(fc-^/l=F) a: ^+-^ dk 

^ ++ {k)+^ ++ {\Fv^ k 2 



k 



+ e 



k 



-i{k+^l-k 2 )x 



+ e 



i(k— vi— k 2 )x 





k 2 


) k 




> y/l- 


k 2 


k 





dk 



where we split I into 
interval. Similarly 



U 



75' 1 



ip(x = y ) 



V2 J 



+ e 



J(k+\/l-k 2 )x 



VT-^k 2 . 

and change variables in the second 



-i(k-VT-k?)x 



-i(k+VT-^)x 



+ e 



i(k—\/l—k 2 )x 





k 2 


) " 




) Vi- 


k 2 


k 





Vi-^k 2 



dk. 



Since for k € ij, —k + \f\ - k 2 £ [0, 1] and k + y/l - k 2 6 [1, \/2], each 
of the four exponentials in the above integrals takes arguments in disjoint 
intervals. As a result continuity is equivalent to the equality of the 
coefficients of the exponentials, i.e. 



tj} (7T (k)y/ 1 - k 2 + (yi - A; 2 ) k = 4> aT {k)^l-k 2 + tp Ta [yi-k 2 ^j k 
for k € 

For the boundary condition ([5]) we again use ([6]) to get 



(14) 



7 i(k+VT^k 2 )x 



k-y/l-k' 



+ e 



-i(k-Vl-k 2 )x 



k+^l-k 2 



+ e 



-i(k+Vl-k 2 )x 



k-^l-k 2 



+ e 



i(k-VT-k I )x 



k + Vl-k 2 



- 4> ++ + ^ ++ 

tp — tp" 



+ 2c^ ++ 
- 2c^~ + 
+ 2c0" 
- 2cijj + 



dk 
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on the diagonal quadrants. Splitting I into 
following equalities of transforms 



V2 



u 



we have the 







+ 2c_ 
+ 2c_ 

-i>- 

-2c + 

- 2c+ 
-ic 



+ 4> 



++ 



(k)^i-k 2 + 



++ 



4> ++ (k)y/l- k 2 + 4> ++ U/l - k 2 ) k 
(k)y/l -k 2 + -$—+<$-- 
4>—(k)y/l - k 2 + U/l - k 2 ^ k 



+ V" 



{k)vl-k 2 + 



4>- + (k)Vi-k 2 + 4> + - U/l - k 2 ) k 

4r + + tp~ + 



(kWi-k 2 + 



(yi - ^ k 

(Vl - A; 2 ) k 

(yi - k 2 ^ k 
(yi - k 2 ^ k 



(15) 



(16) 



(17) 



4> + -(k)y/l- k 2 + ip- + (y/l - k 2 ) k 



(18) 



where c± = —ic [k db \/l — k 2 j and again k € As we need to con- 

sider transforms evaluated at k and y/l — k 2 these equations are certainly 
not linear in the variables defined in however, extending the definition 
of £ and % to 



/ ^± + (k)VT^ \ 

^ + (y^w) k 
\ (yi^w) k ) 



Xij 



I ^±-(k)vryw \ 

$tr (v 7 !^) k 



(19) 



k G we see that 



with 



can be written as 



M 



N 







Mia , xu 


= N X u 




1 


l + c_ 


C- 





° ^ 


— c_ 


l-c_ 
















1 - c_ 


— c_ 


\ 








c_ 


l + c_ / 


/ 


l + c 4 








c+ \ 







1 + C+ 


c+ 










-c+ 


1-C+ 





V 


-c+ 








1-C+ / 



(20) 



Using this new definition of £, x the solutions of the linear system (|13|) 
are simply doubled — r becomes the permutation (13) (24) on C 4 and the 
eigenspaces of r are now 2-dimensional. However, since k e yi in (H2J) 
we still get the same set of solutions: the map (J7|) should be considered as 
a map from L 2 ( J | / ) ® (V © V) where the first copy of V is spanned by 
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elements of the form F % {k\,k%) and the second copy of V is spanned by 
elements of the form F % (&2,&i). In this way we see that (|20p is a set of 
linear equations on V © V with the momenta swopped in the two copies of 
V . But the solution of such a problem is exactly what is described in §4 of 
[2J. Consequently, W describes the complete solution space of (fTTfl |2"U|) and 
(JBJ) describes all the eigensolutions of H c with real momenta. 



3.2 Explicit expressions for eigensolutions of H c 

We begin by recalling the definitions of the one and two particle solutions 
described in [2J; the definitions are slightly changed for simplification. 
The scattering wave solutions {ip l }^ =1 are defined in exactly the same way 
as in [2]. We also need the solutions 

where now j G {1, . . . , n} and ip n+1 = ip 1 (ie. the <\? are not linearly in- 
dependent; however, this allows us to define solutions 'mod n'). In order 
to be consistent with the above notation we need new notation for the one 
particle solution with a discontinuity on the diagonal 

ti x . u) _ J s'm(kxi) : x { G Qa, Xi > yi or a?* G i / j 

\ (1 - n) sin(A;xi) : ac, G Qs, X* < y» 

(with the analogous definition for the second particle (y)). The two particle 
states are then 

U% = ^(x,k s ).^(y,k t )\ 

I. J ije{0,...,n} 

M=^(x,*; a )-^'(i/,*t)i f , 

*. J i jG{l,...,n} 

= k a ) ■ av, h) - t(x, k s ) ■ <p( y , fct)} ie{1> ... in} , 

with (a,t)€ {(1,2), (2,1)}. 

In the following we describe the solutions G % spanning W splitting them into 
symmetric and antisymmetric solutions. Note: here symmetry, antisymme- 
try, is with respect to the exchange of particles (x,y) «-> (y,x). 

1. The n 2 antisymmetric solutions 

{*S-*S}.. Fft 

*. J i,j6{l,...,n.} 

2. The symmetric solutions 

(a) with support outside the diagonal 

+ \ 

I 12 21 /j,j6{l,...,n},|i-i| mod n>2 

of which there are n 2 — 3n, and 
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(b) with support on the diagonal 

U 2 - %1 + ^1 (*« + *£) - ^ (*« + *$) ) 

I c c J ie{l,...,n} 

of which there are n. 

Here \i — j\ mod n > 2 is to be understood as: the distance between i and 
j, as elements of Z/nZ, is > 2. 

The 2n 2 — 2n solutions described above correspond to the solutions described 
in §4 of [2] in the following way: the n 2 antisymmetric solutions correspond 
to the antisymmetrisation of the solutions from §4. 1 along with the solutions 
from §4.2 of [2]; the re 2 — 3n symmetric solutions with support outside the 
diagonal correspond to the symmetrisation of the solutions from §4.1.1-2 of 
[2]; and the n symmetric solutions with support on the diagonal correspond 
to the symmetrisation of the solutions from §4.1.3 and the solutions from 
§4.3, of [2]. 



4 Conclusion 

As stated above, one of the motivations for studying this problem is that 
the equivalent problem on the line is completely solvable for an arbitrary 
number of particles. It is of interest to know how the geometry of a star 
graph effects this solvability. We see that for two particles the (real mo- 
mentum) eigensolutions are, somewhat analogous to the line, all essentially 
described by finite sums of products of one particle solutions (on the con- 
figuration space cut along the diagonal). Clearly it would be interesting to 
consider three or more particles and observe whether the crucial property 
of the Fourier transforms, viz. the boundary conditions amount to a linear 
system on some finite dimensional vector space, is satisfied. 
Another interesting problem is to investigate whether there exist any solu- 
tions with non real momenta. We do not attempt a rigorous examination 
here, rather we merely mention some relevant points. For two particles on 
the line with an attractive potential (c < 0) there appears a branch of the 

continuous spectrum 



, oo j of finite multiplicity [T] . For this to appear 
in the case of the star graph with n > 3 it should be associated with solu- 
tions which have discontinuous derivative on the diagonal, i.e. the solutions 
described in §4.3 of [2] 

G* = *| 2 - ¥ 21 + Vh. ($«* + d>|°) _ (*« + $f 2 ) . 

It can be shown that on the diagonal quadrants this solution has the form 
G 1 \q.. = n(smk'\xi-yi\smk(xi + yi) -sin k\xi-yi\ sink' (xi + yi) - 

— cos k\xi - yA sink'(xi + yA + — cos k'lxi - yA s'mk(xi + yi) ) 
c c J 

with k = (k\ + k<2)/2 and k' = {k\ — k<i)j2 momenta along the directions 
Xi — Vi and x\ + y%. In the case of an attractive potential (c < 0) the 
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finite multiplicity branches of continuous spectrum appear when one of these 
momenta, say k, is set to k = ^ so that the solution becomes exponentially 
decreasing along the direction X{ — y%. Varying the other momentum, as 
a real parameter, we get a branch of the continuous spectrum starting at 
(y) < 0. In the case of G % this gives 

Gi \n-k=i° = in(-e^ Xt ~ y ^s'mk\xi+yi) + 



\2h' 1 

sinh§(xj + y ii 



2k' 

cos k'lxi — yA + sin k\xi — yA 

c 



with the first term exponentially decreasing. However, the second term is 
exponentially increasing with no apparent (to the author) way of removing 
it. For this reason it appears probable that such branches of continuous 
spectrum do not appear for the n > 3 star graph. A related question is 
the description of the deficiency elements for the symmetric operator H c= q 
restricted to functions vanishing along the diagonal D, It would appear 
that these are related to the basic solutions described £[2] (with eigenvalue off 
the spectrum). Such an operator theoretic description of the problem would 
be highly desirable and is indeed one of the motivations for this investigation. 
Such a connection would allow for the application of standard techniques 
from extension theory to describe eg. the greens function for this problem. 
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Appendix 



Here we define {ej}" =1 to be the basis of C n in which 



P 



1 

n 



(I 1 
1 1 

V = 



and {fi}2=i to be the basis of C n in which 

S = 2P - I = 



1 



where the lower right component of the above matrix is in M n _i(C). Unless 
otherwise stated we use the basis {fi} in this section. 

Since the transposition r commutes with S and ir we may split f)13|) into the 
±1 eigenspaces of r. In this way we see that describing all solutions of (fT3|) 
amounts to finding the kernels of 



p±(x,x) = K x (xS±Sx) ,7T 1 ix-x) 
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on the Hilbert space M„(C)©M n (C)-we consider M n (C) instead of M n (C 2 ) 
since the eigenspaces of r are one dimensional. Defining 

Q±(x,x) = (xS±Sx, x -x), 
^(x,x) = (W,^*) 

along with the subspace 

K± = {x e ker x (Q±) ; Q±x £ ker (iT 1 ) } 

it is clear that we have the decomposition 

ker(P ± ) =ker(Q ± )©K ± . 

We see immediately (using {fi}) that 

™ - ((°, ?U! ")) m 



_b o y ' v b o , , 

^« + ) - ((; -/)) <-) 

ker(Q-) - lY" ° W" ° Y) (23) 



In particular null (Q+) = 2(n — 1) and null (Q-) = (n — l) 2 + 1. 

We describe K± by first finding ker (n- 1 -) n ran(Q-i-) and then finding the 

inverse image with respect to Q± in ker -1- (Q±). For this we need 



We first describe -fT + . By definition 

ker (U 1 -) = {{a ® ei)? =1 , <e< ® e;)" = i) 



in particular ker (n- 1 ) consists of symmetric matrices so 
ker(n^)nran(Q + )c((^ ° ) , ( »' ° 

Consequently, we seek matrices diagonal in the {e^} basis which commute 
with P; but the only such matrices are the scalar matrices so 



ker (u^ n ran (Q+) = (2d, 2a!l) . 
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The preimage of this in ker -1- (Q+), from (I22p . is 

a + a' \ f a-a' 

-(a-o')I / ' V 

so that dim (K+) = 2. 
For K_ we have 

ker(n ± )nran(Q_)c (Y £ * V 





-(o + o')I 



a b' 
b' -A 



again using the fact that ker (n- 1 -) consists of symmetric matrices. Here the 
repeated b (&') are understood to be transposes of one another. We should 
seek solutions of the form 



i 

so that 



Cj &i 09 6j 



a b 
b A 



es «> e.; 



a y 

b' -A 



X = X 



+ X =J2 °i 

i 

with c'l = Ci + d i . By construction 
= P ± X "P ± 



e,- » e,- 



2a 6 + 6' 
6 + 6' 



e,: <K> e,- 



4(e«)®(i> 

-£(e*We<* 

^ ^ e< (8 ei - - e i ® e i c " + c " " - E c 'k 



This implies that c" = q + c' { = so that 

b 



ker n 1 nran(Q_) C 



b A 



-b 
-b -A 



(27) 



which in turn gives 



P X P 

^E cfc E 



d ej 



h3 



or Y^k c k = 0- in this way we see that ker (n J -)nran (Q-) is n— 1 dimensional 
and is spanned by 



^ Ci ei ® e,, - ^ Q ej <g> e» J , ^ c fc = . 

v i i / k 



(28) 
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We describe the inverse image: using (|24h we see that the inverse image of 
([27D in ker- 1 (Q_) is 

K ~ C (( A b ,2 )•{-„ J/2 )) ' 
Suppose we have solutions of the form (|28p then we claim that 

££«*®E<* = (2 o) 

This is seen to be true since 

(&) = (° l) /l= (^ Cie ^ ei ) /l = ^^ Qei ' 
Consequently, the inverse image of Q~ in ker -1- (Q-) is spanned by 

I E c * e * ® e * + - E ei E C J e i + ~ E c * ei ® E e J ' 

\ i j i i j 

- (H ei ® ei + ^ E e i ® E C J e i + ~ E Ci ei ® E e J 

i i j i j / 

with Cfc = 0. In particular it also has dimension n — 1. 

In summary, the solutions of (fT3|) split into four subspaces: ker (Q_), ker (Q+), 

K_ and of dimensions (n — l) 2 + 1, 2(n — 1), n — 1 and 2 respectively. 
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